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1. Introduction 

In order to test the theory of strong interactions (QCD) it has long been common 
practice to confront experimental data with perturbative predictions and then extract 
fundamental parameters such as the strong coupling a s . A consistent extraction of 
a s from different experiments and different observables points at the correctness of 
assuming QCD as the underlying theory. 

Since QCD is a theory where non-perturbative effects play a clear role one has to 
first question the feasibilty of such an approach and review the role of perturbative QCD 
(pQCD). The fact that one can make perturbative predictions to any extent is indeed a 
blessing and in the case of many observables factorisation theorems allow the separation 
of pQCD from non-perturbative effects. Simply stated a factorisation theorem implies 
for an observable V 



In the above we indicate by x any set of phase space variables on which the observable 
may depend and by Q 2 the hard momentum scale. O is typically a universal operator 
matrix element which contains non-perturbative information and whose value can be 
extracted from experiment and tabulated, for example a parton distribution measure. 
C\ and C 2 are perturbative coefficients calculated by using Feynman graphs to represent 
the hard process. 

The dots indicate the limitations of the perturbative approach. The dots in 
the brackets tell us that there are missing terms in the perturbative series since any 
perturbative predictions (even resummed ones) are missing higher order pieces which 
cannot be easily computed. While this is just a technological limitation in the sense that 
higher order calculations become rapidly prohibitive computationally, it is a common 
misconception that given adequate computing power one could have a well defined 
perturbative answer that would be arbitrarily accurate. 

That this belief is ill-founded can most easily be seen by examining a particular 
type of contribution (renormalon graphs) (see for examples and further references |l|) 
which give that for large enough n the coefficients C n diverge factorially. This does 
not mean that the concept of a perturbative expansion is useless. What it implies is 
that pQCD gives an asymptotic approximation to the observable V. The error made in 
truncating the perturbative expansion is of the order of the first omitted term. Therefore 
one should carry out the perturbative expansion till the smallest term is reached. The 
optimal error which is of the order of this smallest term is thus an unvoidable inherent 
ambiguity in the perturbative expansion and one finds 



where A and p depend on the observable and 5V represents the uncertainty of the 
perturbative estimate . The coefficient A may contain dependence on one or more phase 
space variables on which V depends. 




(1) 
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Also in Eq. (1) there are dots to indicate that the simple factorisation may not 
be correct at some level, e.g higher twist. These contributions are also typically power 
behaved and the expectation is that in order to obtain an unambigous answer the 
ambiguity of the perturbative (PT) part is cancelled by an ambiguity in the definition 
of the non-perturbative (NP) higher twist operators. Assuming further that the size 
of the ambiguity is proportional to value of NP terms (condensates) one can use the 
ambiguity in the PT piece to model NP effects which are of a power law type. In other 
words where there is a large power behaved ambiguity in the PT series it can be treated 
as a signal that NP effects are damaging in proportion. This assumption is known as 
ultraviolet dominance of higher twist. 

Having arrived at the notion that the ambiguity of the perturbative expansion is 
relevant to the study of power corrections, one is still faced with the task of computing 
it and comapring with experimental data. This is conveniently done in the framework 
of the dispersive approach we shall next discuss. 

2. Dispersive approach 

This method was formulated by Dokshitzer, Marchesini and Webber |2| and relies on 
the assumption that the power correction to an observable is given by evaluating the 
perturbative renormalon graphs and estimating the ensuing ambiguity of the PT series. 
To investigate this method one observes that inserting renormalon bubbles on gluon lines 
in a lowest order Feynman graph for a given observable V results in the appearance of a 
running coupling constant. If this gluon is an internal line inserting renormalon bubbles 
gives for an observable V 



with k 2 representing the gluon virtuality. The role of the renormalon bubbles is simply to 
give rise to the running coupling in the integrand above. Strictly within a perturbative 
context the above integral gives a divergent result due to the divergence of a s at the 
Landau pole. This is just the factorial renormalon behaviour discussed previously. 

To generate power corrections from the above contribution the assumption one 
makes is that the QCD coupling a s instead of diverging is a well behaved universal 
function down to small scales, i.e infrared finite. Further knowledge of the exact 
behaviour of the coupling is not required for us to proceed. This is certainly not an 
assumption that can be justified from first principles. However it leads to simple results 
that can be easily tested against a wealth of experimental data. This exercise has been 
carried out for different observables and so far it seems that experimental data are in 
support of this notion of a universal IR finite coupling. 

Assuming that the coupling is analytic except a branch cut across the negative 
real axis which is a requirement of causality, it follows that one can write the formal 
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dispersion relation 

r°° dm 2 

<k 2 ) = - -^uPsim 2 ) (4) 

where the spectral function p s represents the discontinuity across the cut, 

Ps(m 2 ) = -^Disc {a s (-m 2 )} = -L {a s (m 2 e™) - a s (mV i7r ) } . (5) 

Non-perturbative effects at long distances are expected to give rise to a modification 
in the strong coupling at low scales, Sa s , which generates a corresponding modification 
in the spectral function : 

5p s (m 2 ) = -'-Disc {5a s (-m 2 )} . (6) 

From Eqs. ([3|-^]) it follows that the non-perturbative contribution to V, trigerred by 
a virtual gluon with renormalon insertions and an NP modified coupling can be written 

as 

f°° dm 2 

SV mrt = / 5 Ps (m 2 )^[V(m 2 ) - V(0)] (7) 
Jo m 

A similar equation holds for the real emission contribution and hence in what follws we 
shall remove the superscript and talk about the sum of real and virtual pieces. 

The quantity V(m 2 ) represents the 0(a s ) virtual calculation performed with a finite 
gluon mass m. The gluon mass is just the dispersive variable in Eq. [| and does not 
imply that one works with a massive gluon field and spoils gauge invariance. 

After some further manipulation one gets the following non-perturbative 
contribution to V 

sv= rM 5as{rn 2 )g{m 2 /Q 2 ) (8) 



m 

where, setting m 2 /Q 2 = e, 

e(e) = -^rDisc{V(-e)} . (9) 

Since 5a s (m 2 ) is limited to low values of m 2 , the asymptotic behaviour of 5V at large 
Q 2 is controlled by the behaviour of V(e) as e — *• 0. We see from Eq. |9| that no terms 
analytic at e = can contribute to 5V. On the other hand for a square-root behaviour 
at small e, 

V-ay^Ve *V-=--#> (10) 

2tt 7T Q 

where we defined the coupling moment 

r°° dm 2 

Ax = / — 5- m 5a s (m 2 ) (11) 
Jo m 

We can express Ax in terms of the average value of a s in the infrared region, as follows. 
We substitute for 5a s in Eq. (pi]) 

5a s (m 2 ) ~ a s (m 2 ) - a p s T (m 2 ) , (12) 
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where a^ T represents the expression for a s corresponding to the part already included 
in the perturbative prediction. If the perturbative calculation is carried out to second 
order in the MS renormalization scheme, with renormalization scale ft\, then we have 

af(m 2 ) = a s {ft\) + [Mn(^/m 2 ) + k] a 2 s {ft\) (13) 
where for Nf active flavours (Ca = 3) 

= UCA-Wf = (67-3^)^-10^ 

12vr ' 36tt 1 j 

The constant k comes from a change of scheme from MS to the more physical scheme 
in which a s is preferably defined at low scales. Then above some infrared matching 
scale ft we assume that a s (m 2 ) and a^ T (m 2 ) approximately coincide, so that 

Ai ~ — / dm (a s {m 2 ) — a s (ft\) — [b ln(/i 2 /m 2 ) + k] a 2 (/x 2 )) 



7T 



— ft (ato(ft) - a s {ft\) - [bln^l/ft 2 ) + k + 2b] a 2 s (ft\)) , (15) 



7T 



where 



1 



rfti 

ao(ft) = — / a s (m 2 ) dm . (16) 
/^i Jo 

Studies of event shapes in e + e~ annihilation suggest that «o — 0.5 for ft = 2 GeV, 
which translates into a value of A\ — 0.2 GeV for Q ~ 20 — 100 GeV. 

3. DIS event shapes 

Here we shall concentrate on the phenomenology of power corrections to event shape 
variables measured by the HI and ZEUS collaborations. This study is carried out in 
the Breit frame and the current hemisphere H c is used to define the variables, freeing 
us from contamination by the proton remnant jet. 

Some of the observables studied are defined below 

• Current jet thrust This is constructed by summing the modulii of hadron 
components along the photon axis in the Breit frame. Two normalisations are 
possible, to Q/2 on one hand and to the energy in the current hemisphere (provided 
this exceeds a certain preset minimum). 

r Q = l-T Q = l-2 Eaeg ° |Pa -" 1 (17) 



r E = 1 - T E 



Q 

Current hemisphere jet mass 

2 



p= /Q 2 

\a£H c J 



(19) 
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• C-parameter The C-parameter is 

C = 6 Yl \Pa\\Pb\sm 2 e ab /Q 2 . (20) 

a,beH c 

There are still other possible definitions and variables studied in addition to the 
ones mentioned above (an important one being the jet broadening) but for illustrative 
purposes it is enough to restrict the discussion to the variables mentioned. 

4. Power corrections 

4-1. Naive results 

The power corrections for different event shape variables are given by following the 
dispersive recepie. One first computes the one-loop graphs with a finite gluon mass eQ 2 
For all the variables mentioned above the leading non-analytic behaviour at small e is a 
\fl singularity. 

This leads to 1/Q power corrections of the form of Eq. |10| 
The coefficients a v differ from variable to variable and one gets 

a p = —2, a T = —4, a c = — 67r (21) 

The constant A\ is expected to be approximately universal (to about 20%)as observed 
in e + e~ data and with similar values to those found there. 

One point that needs to be made in the DIS context is that the leading power 
corrections to most event shapes are independent of the Bjorken variable x. This comes 
about since the 1/Q effects in question are trigered by a soft gluon, the emission of 
which does not change the momentum fraction carried by the incoming projectile. 
An exception to this is the jet broadening variable where an interplay between this 
emitted soft gluon and x dependent perturbative radiation gives rise to an x dependent 
power correction. A complete estimate for the broadening power correction relies on a 
resummed perturbative calculation which shall be published shortly f|J. 

Another point that needed to be addressed was one first raised by Nason and 
Seymour ||. They pointed out that gluon decay can only be mimimcked by a massive 
gluon provided one can integrate inclusively over the gluon offspring. In other words 
the machinary of the dispersive approach should only strictly be applied to inclusive 
observables which do not register offspring gluons. Event shapes are not such variables 
as the offspring gluon momenta affect the event shape. Hence one needs an improved 
approach where one explicitly treats gluon decay properly. 

4-2. Two-loop improvement 

The two-loop calculations that take gluon decay into proper account were first performed 
for e + e _ observables |j) Following this similar work was performed for DIS observables 
0. The conclusions could at first sight be described as surprising. Taking two loop 
effects into account meant a simple rescaling of the old one-loop coefficients by a 
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Figure 1. Summary of HI results for fits to event shape means. Figure taken from 
Ref. 0. 



universal factor A4 = 1.49 for rif = 3. Hence the ratios of power corrections to different 
observables are still given by the ratios of the one- loop coefhcents ay. 
This surprising result came about clS db result of 

• Universality of soft gluon radiation. 

• Common geometrical properties of event shapes -they are all linear in the emitted 
transverse momenta. 

4-3. Differential distributions 

Another area of investigation at HERA is the study of event shape differential 
distributions and power corrections. Here while data has been available for a while 
there was a lack of theoretical estimates. Although there are Monte Carlo programs 
to make fixed order estimates at 0(a 2 s ) this is not sufficient in the above case. The 
presence of large logarithms as big as a™ ln 2n_1 y at n th perturbative order means that 
a fixed order expansion is not very useful in the small V region. 

As in e + e~ one needs to resum all orders till next-to-leading logarithmic accuracy 
to describe the data. This resummation program has been started |IJ and results for 



many variables are soon to be published ||] . However a preliminary comparison to data 
has already begun. 

The power corrections to different event shape distributions result in simple shifts of 
the perturbative spectra, the amount of the shift is exactly equal to the 1/Q correction 
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• Stat-only fit, Q>9 GeV 

* Stat + Sys fit, Q>9 GeV 
(x 2 = Xmin + 4 contours) 
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Figure 2. ZEUS contour plots for event shape means. Figure taken from || 



to the corresponding mean values. 

where a v for the different variables are quoted. 

Again the jet broadening is different and does not show a simple shift of the PT spectrum 
giving instead a variable dependent shift as in e + e~. 

5. Fit results 

The experimental results from HI and ZEUS (see Figs. 1 and 2 ) have been available for 
quite some time for mean values whereas fits to distributions are on the way. For the 
mean values the fits are best summarised by plotting error ellipses in the a s ,ao plane 
where do is the parameter introduced in Eq. [15|. 

Some brief comments about the results are in order. In the HI results one should 
note that the jet mass ellipse is shifted and falls in line with expectations from other 
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variables provided one neglects hadron masses (the relevant ellipse is labelled p ) m 
keeping with the theoretical approximations made. For a fuller discussion of this point 
see [11]. In the ZEUS plot one notices the effect of large systematic errors. ZEUS have 
also studied the x dependence of power corrections and as noted before in this article 
the power correction to the jet broadening variable displays a clear x dependence which 
comes from the anomalous dimension matrix J3J . 

In addition one remarks that that the values of the oiq parameter obtained from 
DIS event shapes has similar values (0.4-0.5) as for e + e~ variables and this seems to 
support the hypothesis of a universal infrared finite coupling. The spread in a s values 
however is still uncomfortably large and is a potential area for improvement. 

Finally turning to differential distributions one can use a nominal value of a s = 
0.118 and «o = 0.5 and the resummation of [0. The curves in Fig. 3 are produced by 
binning like the HI data which is responsible for their non-smooth nature. Preliminary 
indications are positive as shown by Fig. 3 but detailed fits are awaited. 
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